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We consider three-body recombination into deep dimers in a mass-imbalanced two-component 
atomic gas. We use an optical model where a phenomenological imaginary potential is added to the 
lowest adiabatic hyper-spherical potential. The consequent imaginary part of the energy eigenvalue 
corresponds to the decay rate or recombination probability of the three-body system. The method 
is formulated in details and the relevant qualitative features are discussed as functions of scattering 
lengths and masses. We use zero-range model in analyses of recent recombination data. The 
dominating scattering length is usually related to the non-equal two-body systems. We account 
for temperature smearing which tends to wipe out the higher-lying Efimov peaks. The range and 
the strength of the imaginary potential determine positions and shapes of the Efimov peaks as well 
as the absolute value of the recombination rate. The Efimov scaling between recombination peaks 
is calculated and shown to depend on both scattering lengths. Recombination is predicted to be 
largest for heavy-heavy-light systems. Universal properties of the optical parameters are indicated. 

We compare to available experiments and find in general very satisfactory agreement. 


I. INTRODUCTION 

The Efimov effect is very well established and the en¬ 
ergies and radii of the infinitely many three-body bound 
Efimov states occur in geometric series corresponding to 
well defined scaling factors depending on the masses of 
the particles. These states have been searched for and 
found, but only indirectly as enhanced decay probabil¬ 
ities when the two-body interactions are tuned to spe¬ 
cific values. Most of these findings are for three identical 
atoms where the geometric scaling factor is 22.7 inde¬ 
pendent of the atomic mass, see for example the exper¬ 
iments i-El- The occurrence conditions for the Efi¬ 
mov effect are that at least two of the three two-body 
subsystems simultaneously interact with infinitely large 
scattering lengths, or equivalently with bound states of 
zero energy. 

The necessary fine-tuning is most easily achieved when 
at least two of the particles are identical. Then obvi¬ 
ously at least two subsystems are identical and therefore 
simultaneously fine-tuned in an experiment. If further¬ 
more one of the particles is distinctly different, the scal¬ 
ing factors can differ enormously for the resulting Efimov 
spectra. When the system has one light and two heavy 
particles, the scaling factor, or the energy spacing, is very 
small. This is a huge advantage since many energies then 
are within the experimentally accessible range. The op¬ 
posite case of one heavy and two light particles has cor¬ 
respondingly large energy spacing and fewer accessible 
states. 

The two-component gases, where mass-imbalanced 
three-body systems are likely to recombine in the decay 
process into a deeply bound dimer and a third parti¬ 
cle, are more difficult to prepare in sufficiently cold tun¬ 
able states. However, several systems are now investi¬ 
gated experimentally [20l - l^ and probably more are in 
the pipeline. 

It is therefore very timely to investigate these processes 


theoretically [2^ in more details and more systematically 
than done before, where one of the foci has been on 
the special case of vanishing scattering length between 
the two similar particles, see for example [^. We shall 
in this paper only consider negative scattering lengths 
where the final state in the recombination product is 
a deeply bound dimer, unavailable in usual zero-range 
three-body models. The theoretical formulation is sim¬ 
plest for three equal masses, since then the single, and 
necessarily large, scattering length is the only parameter 
which therefore completely determines the properties of 
the process. The two-component system is much more 
complicated, since two scattering lengths and one mass 
ratio are the necessary parameters to characterize the 
system. 

The scaling factor is very well known as function 
of masses in the limit of very large scattering lengths 
[2^ [2^ , that is when both are infinitely large and when 
only that of the non-equal particles is large while the 
third between identical particles is vanishingly small. 
These two extreme limits of different scalings between 
energies would lead to different periodicities in the en¬ 
hanced recombination probability as function of the scat¬ 
tering lengths. Any intermediate energy scaling is also 
possible, and a smooth transition between these limits 
arises by varying the third scattering length between zero 
and infinity. 

One established method is based on a multiple scat¬ 
tering approach where the final state somehow is simu¬ 
lated although never directly present. Recently a novel 
suggestion appeared to use an optical model to describe 
absorption as in nuclear physics and for light waves in 
materials. This formulation can only provide decay, or 
absorption, probability without any details of the final 
states populated in these processes. However, this is all 
we need to describe the measurable recombination prob¬ 
abilities. In the present report we shall elaborate and 
explore the properties of this optical model. 
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The overall purpose of the present paper is to discuss 
three-body recombination processes in atomic gases. The 
emphasis is on two-component atomic gases where mass- 
imbalanced recombination processes are prominent and 
probably dominating. This shall be done by introducing 
the details of a simple optical model where the real part 
is the lowest adiabatic hyper-spherical potential, and the 
imaginary part is a square well in hyper-radius. 

The initial zero-range potential must be regularized to 
provide physical results. This is achieved by replacing the 
diverging real part with a constant at small hyper-radii 
where the imaginary part also is finite. This constant 
is unimportant for the calculated recombination rates. 
Such an optical model formulation is phenomenological 
by nature, but we shall indicate a possible relation to 
universally determined values of the model parameters. 

The paper is structured with basic definitions in 
section II, while the focus on computation of mass- 
imbalanced three-body recombinations is formulated in 
details in section III. The choice and dependence on phys¬ 
ical as well as model parameters are discussed in section 
IV. Precise comparison to available experimental results 
is the content of section V. Finally, we summarize a num¬ 
ber of conclusions in section VI, where we also suggest 
interesting future research projects. 


II. NOTATION AND BASIC INGREDIENTS 

In this section we sketch the details necessary to cal¬ 
culate the potential which provides the wave functions 
for recombination computations. The first part describes 
the adiabatic expansion method and the dependence on 
masses and scattering length parameters for a short- 
range real potential. The second part discusses the ex¬ 
tension to an optical potential by addition of a complex 
term. 


A. Hyper-spherical expansion 


We utilize the formalism developed in to 

treat the 3-body problem at low energies. The parti¬ 
cle coordinates and masses are and m^, respectively 
for i = 1,2,3. We describe the three-body system 
by hyper-spherical coordinates where the all-important 
hyper-radius is defined by 


2 ^ 2 , 2 J2^<k'^irnk{r^-rk)^ 
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with an arbitrary normalization mass, m, which in the 
present work is chosen to be the mass of the lightest par¬ 
ticle, unless something else is explicitly stated. The first 
term in the hyper-spherical expansion is known to give 
fairly accurate results [2^. This is provided the exci¬ 
tation energy is less than the energy difference between 
states build on first and second adiabatic potential. If 
we extract the radial phase-space factor, the total wave 
function, 'k, is approximated by 

«'(p, flp) Ri /(p)$(p. Op), (4) 


where Op denotes the other five hyper-spherical coordi¬ 
nates. The hyper-radial wave-equation to determine / is 
then given by [2^ 



2mE \ 
^ ) 


f{p) = 0 , 
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where z/ is a function of p. The procedure to determine 
V is described in details in (2^. Briefly, the adiabatic 
hyper-spherical expansion method amounts to fixing the 
hyperradius and solving the three Faddeev equations in 
the space of the five dimensionless coordinates. Op. All 
partial angular momenta are zero corresponding to s- 
wave spherical harmonics solutions for the angles de¬ 
scribing the directions of x^ and y^ in all three Jacobi- 
coordinate sets. We are therefore left with the three cou¬ 
pled Faddeev equations describing the relative sizes of Xi 
and yi- For zero-range two-body interactions these differ¬ 
ential equations can be reduced to three linear algebraic 
equations with non-trivial solutions only when the cor- 
resp onding 3x3 matrix, M, has vanishing determinant 

m 


det M = 0, 


( 6 ) 


which determines the angular eigenvalue, v. The matrix 
elements of M are given by [^ . 
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The scattering length, ajk, is for the interaction between 
particles j and fc, and the angles, (pij, are given by 




i = arctan 
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where {i,j, k} is a permutation of {1, 2,3} 

For the simplest case of 3 identical bosons, Eq.dH) re¬ 
duces to (Mil — Mi 2 )(Mii -I- 2M12) = 0, where we only 
have to distinguish between diagonal and non-diagonal 
terms. The second factor equal to zero is equivalent to 
the well-known equation for v{p) 
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FIG. 1: Visual illustration of the numbering convention used 
to calculate the adiabatic potentials. Particle 2 is distinctly 
different from the identical particles 1 and 3. 


where all indices are omitted for this case of identical 
bosons. The other case of interest is two identical and 
one distinguishable particle. We label as illustrated in 
fig. n that is particles 1 and 3 are identical (toi = m 3 , 
0-12 = 0,23, Mi 2 = M 23 , and Mn = M 33 , M 13 ^ Mu). 
The characteristic equation then reduces to 

(Mil ~ Mi 3 )(MiiM 22 — 2Mi2 + M22M13) = 0 , (11) 


where the last factor provides the general physical solu¬ 
tion, which receives contributions from all three interact¬ 
ing pairs. Non-interacting particles 1 and 3 correspond 
to vanishing scattering length, 013 = 0 , where only two 
interactions contribute. Then M 22 —>■ 00 and the sec¬ 
ond factor in Eq. m reduces to Mu -|- M 13 = 0 , which 
explicitly amounts to 


i^cos 



2 sin(t/[(/)i 3 - f]) 
sin(2(/)i3) 


y/jHai2 


■ sin ^ 


z/7r\ 

T / ■ 


( 12 ) 

In all the calculations of v{p) we only use the small¬ 
est non-trivial solution to the Eqs. (nni), di) and (Ha. 
This selects the lowest adiabatic potential in the hyper- 
spherical expansion. The results obtained for v depends 
on the ratios of hyper-radius to the different scattering 
lengths. 

When Oi —> 00 the purely imaginary solutions, v = 
vq, become independent of p, where t'p depends on the 
masses and the number of contributing interactions. The 
radial equation in Eq. ([S]) has then infinitely many bound 
solutions, the Efimov three-body states, with energies, 
En, related by 


^”+1 _ „-27r/|i/o| = J_ 

En 


(13) 


where we generally refer to s as the Efimov scaling. We 
emphasize that s in particular depends on the number of 
contributing subsystems. In practical comparison with 
measurements this means that we can expect s to take 
any value between the values of the two extreme limits 
of two or three contributing subsystems. 


B. Optical potential 

We are interested in modelling the recombination for 
negative values of all the scattering lengths, that is in 


the regime where no two-body bound states exist within 
the present zero-range model. We therefore need to in¬ 
troduce the final dimer states populated in the final state 
of the three-body recombination process. However, the 
only information we need is the rate of population, or 
equivalently the rate at which the three-body system dis¬ 
appears. No details of the final states are required, and 
the optical model is perfect for this purpose [1^. This 
model has recently been employed to describe three-body 
recombination of identical bosons M- 

We add an imaginary part to the adiabatic hyper- 
radial potential in Eq.([5]) leading to non-hermiticity and 
non-conservation of probability. Thus we can describe 
the time-dependent probability reduction as an absorp¬ 
tion process, since recombination is removal of probabil¬ 
ity from the initial three-body system. This is then di¬ 
rectly aiming for calculation of the absorption rate. The 
recombination process occurs when all three particles si¬ 
multaneously are close in space and two can merge into 
a dimer, while the third is necessary to conserve energy 
and momentum. The hyper-radius is an appropriate and 
convenient measure of average distance between the three 
particles, see Eq. ©■ We therefore modify the potential 
in Eq.([5]) at small distances, p < Pcut, by 




-14 


imag 


if P ^ Pent 
if P ^ Pcut 


(14) 


where peut and Vimag are constants characterizing the 
optical potential, V{p). This structure regularizes the 
otherwise diverging, for p —?► 0 , real part of the potential 
by using a p-independent constant for p < peut- 

The solution, /(p), to Eq.© is that of a complex 
square well for p < Peut, that is 

/(p) = Hsin(Kp) , K = ^j 2 m{E - Vipeut))/^^^, 

( 15 ) 

where both k and V are complex quantities. This solu¬ 
tion can conveniently be used as initial condition in the 
numerical solution for an integration starting in peut- 
A schematic illustration of the real part of the poten¬ 
tial is shown for three identical bosons on fig. [21 This 
potential has a short-distance attractive region, a barrier 
of height 0.2877i^/(2TOa^) located at p = 1.46a, and a 
decrease towards zero as 15/(4p^), where the latter can 
be seen from Eq.® with v = 2 obtained from Ea. dTUl) . 
since sin(i 27 r/ 2 ) = 0 in the limit of p —^ 00 . 

The small distance attractive behavior outside Peut, 
but inside the barrier, is proportional to —Ijp^. Such a 
potential produces a number of bound three-body states, 
which only is limited by the finite value of V for p < p^ut ■ 
The energies of these states are related through the scal¬ 
ing in Ea. (n® . Physically we can think of the recombi¬ 
nation process as related to the probability of reaching 
the absorptive small distances by tunneling through the 
barrier. The recombination probability is then substan¬ 
tially enhanced when one of these states has an energy 
equal to the total three-body energy. 
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FIG. 2: A schematic illustration of the real part of the opti¬ 
cal potential for three identical bosons. Both axes use natu¬ 
ral length and energy units of a and /{2ma^). For values 
p > Pent it corresponds to the adiabatic potential for identical 
bosons. At smaller p it has a constant real value correspond¬ 
ing to the value of the adiabatic potential in pent- The blue 
lines correspond to some of the bound states of the —Ifp^ 
potential. The dotted red line illustrates a possible resonant 
state at positive energy. The scale of this illustration is un¬ 
physical, since the value of peut in actual calculations tends to 
be much smaller relative to the value of p where the potential 
has a maximum. 

The scattering length is decisive for potential and 
bound state energies. In particular, for three identical 
bosons we define a^~'^ as the value for which a bound state 
appears with zero energy. The recombination probabil¬ 
ity consequently peaks for the same, very small, energy 
close to zero when a = The parameter peut turns 

out to determine the value a^~\ whereas the strength 
parameter Vimag determines the shape and size of this 
recombination peak as a function of a. 

III. THREE-BODY RECOMBINATION 

The potential described in the previous section now 
has to be used to calculate the rate of probability dis¬ 
appearing in the different three-body channels. We first 
define formally the different channels and rates, then we 
describe how to compute these rates, first by the tra¬ 
ditional method employed for identical particles, and in 
the last subsection we discuss a new method intuitively 
related to the optical model. 

A. Rate equations 

For mass-imbalanced recombination the rate equations 
for the loss-rates are more complicated than for the mass- 
balanced case. For two components with densities ni and 
712, we can have three-body systems with one distinguish¬ 
able and two identical particles. In total we then have 
4 different possible recombination processes, and we can 
in principle find the time derivative of either of the two 
particle densities which leads to 6 different recombination 


coefficients, denoted a, that is 

n'l = — 0221^2’^! I (1®) 

(2) 3 (2) 2 (2) 2 

n2 = -a222?^2 - all2^^2nl - a^2>2«i- (1') 

Here the indices 1 and 2 refer to the two different kinds 
of particles in a two-component gas. The first term in 
Eas. lTBl) and (flTl) corresponds to the recombination co¬ 
efficient for identical particles, whereas the last two terms 
correspond to the two mass-imbalanced recombinations. 

We shall in the following formal derivations often imag¬ 
ine the example of a Cs-Li gas where Cs is particle 1 and 
Li particle 2. Then a^ii 2 corresponds to the recombi¬ 
nation process between two particles of type 1 and one 
particle of type 2, and the rate relates to the change in 
the particle density of type 1. Analogously for other sub- 
and superscripts. 

The rate equations only describe how many particles 
disappear from a three-body system. The final state is 
not specified. We could attempt to further split the mass- 
imbalanced recombinations into two types depending on 
the final structure of which particles form the dimer, that 
is 

1-f 1-f 2^ l-h(12) or 1-hl-f 2^ (ll)-h2. 

However, the optical model does not allow distinction 
between final states, and there is also currently no way 
to distinguish the two channels experimentally. Thus, 
is the full recombination coefficient corresponding 
to the sum of these different dimer productions. 

It is obviously much more difficult to find ni{t) and 
712 (t) from Eas. dTBl) and (IT71) than solving one much sim¬ 
pler equation corresponding to a gas of identical bosons. 
However, we assume the coefficients al*]. are density and 
time-independent, and we are able to calculate these co¬ 
efficients directly from the radial equation in the optical 
model-modified Eq.@. Also experimental analyses are 
more direct by measuring the loss of individual particles 
as function of time. 

We therefore do not attempt to find the full time de¬ 
pendence from Eqs. m and (1171) . It is, however, reassur¬ 
ing to estimate the time-scale at which the numerically 
computed values of afjf, make sense. Let us first assume 
that we have a single-species gas corresponding to the 
equation rii = —The solution is 

^ ‘ ( 18 ) 

\l^ + 2ani^i(^o)(^ - to) 

This square root dependence holds for the single-species 
gas. With typical experimental parameters for the den¬ 
sity [l^, we find that ni is reduced by a factor of 2 over 
a time period varying from about half a second to a few 
nanoseconds as the scattering length changes from -lOOoo 
to -20000ao, where oq is the Bohr radius. We used here 
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the values of obtained in the calculations reported 
in details in a later section of this paper. 

Let us then consider Eq. (fTOl) but with 712 varying slowly 
enough to assume it is constant. If now the 0^22 recom¬ 
bination is dominant, we get the solution 


cients, we need to transform to densities using the iden¬ 
tity Hi = ^. By insertion of this in Eas. (l21l22l) and 
comparison with Eas. (ll6ll7l) this yields the following re¬ 
lation between recombination rates and recombination 
coefficients: 


ni(to) 1 + a[\\ni{to)n2{to){t - to) 

Under these assumptions, we get a half-life ranging from 
a few milliseconds to a few nanoseconds as the 012 scat¬ 
tering length grows from -lOOoo to -lOOOOao- These num¬ 
bers correspond to a Cs-Li gas, where the densities are 
obtained from the experiment [ 2 ^ and the values of 
are from calculations reported later in this paper. 

If the last term is dominating in Eg. (11611 we get anal¬ 
ogously an exponential time dependence, that is 

= exp(-a^ 2 i»^ 2 (^o)(< - to)). (20) 

ni[to) 

Using the densities obtained from the experiment [^ . 
this half-life ranges from about 150 seconds to a few mil¬ 
liseconds as the 012 scattering length increase from about 
-lOOoo to -IOOOOoq. The relatively large half-life here re¬ 
flects the small values of 0 : 221 : which are from calculations 
reported later in this paper. 

We shall return with more discussion on calculated val¬ 
ues and relative sizes of the recombination coefficients in 
Eas. lfTOl) and (fT71) . In general, our results for the rate 
coefficients, ob)., combined with our simple approximate 
solutions give realistic time-scales compared to experi¬ 
mental conditions. When analysing the data in a given 
experiment, a variety of methods are used. Some involve 
a numerical solution of an equation similar to Eg. (1161) 
[^ . while some are more elaborate taking the experi¬ 
mental variation of the temperature with time into ac¬ 
count [ 3 ^ . 

In theoretical calculations we generally find the prob¬ 
ability loss per unit time for a single three-body system, 
denoted Jijk^ which will be refereed to as the recombina¬ 
tion rate. Erom this the probability loss of single parti¬ 
cles involved in the recombination can be found. Let us 
look at the particle loss, corresponding to the terms in 
Ea. (fT51) . The number of particles 1 lost per unit time, 
Ni, from one three-body system is respectively 3J, 2J 
and IJ for the (1-1-1), (1-1-2) and (2-2-1) systems. The 

total number of three-body systems is respectively 

and for the (1-1-1), (1-1-2) and (2-2-1) sys¬ 

tems. Equivalent arguments hold for Eq. (O and thus 
we find the total particle loss per unit time to be 

A = ^ Jii2iV2iV2 - ij22iiViiV| (21) 

N 2 = —^•/222-^2 ~ ~ 2 


®111 = > *^222 = 2'^^^^^222 ( 23 ) 

0^112 “ '^112Ul12 j Q!i12 ~ (24) 

*^221 “ 2 '^22ll221 I Q!221 “ '4221^221 • (25) 

The indices are necessary to distinguish between the four 
different three-body systems. The quantity V^jf.Jijk, is 
then independent of the volume, under the assumption 
that the volume is sufficiently large. The volume used 
in theoretical calculations and the volume corresponding 
to specific experiments are not the same, but 
should be the same and thus it is the quantity for which 
comparisons between experiment and theory is meaning¬ 
ful. The difference between V^f,Jijk and the recombi¬ 
nation coefficients is the factors related to the number 
of three-body systems and particles lost per three-body 
recombination, all of which have been accounted for in 
Eo. (l^51) -Ea. (ES]) . For a non-BEC gas there is also a mul¬ 
tiplicative symmetry-factor corresponding to the number 
of particle permutations Pijk, as described in 13311. b ut for 
a BEC-gas this symmetry factor disappears |33l. l34| . For 
three identical particles Pm = 3!, while for two identical 
particles and one distinguishable particle Puj =2!. 

In the following two sections, two methods of calculat¬ 
ing Jijk for any three-body system is derived. In order to 
ease the notation the subscripts on J and V are omitted. 


B. S-matrix method 

The recombination rate J for a given three-body sys¬ 
tem is defined from the missing probability after scat¬ 
tering on the optical potential in Eo. dHl) . Let us assume 
the three-body wave-function asymptotically is expressed 
in Jacobi coordinates as a three-dimensional plane wave 
normalized in a box of volume, V. We expand this wave- 
function in hyper-harmonic free solutions, that is [ssj] 

■0 = J_gik,cX+ikyy _ 

where kx and ky are the Jacobi momenta characterizing 
the wave function, = {k^ + ky), and denote the five 
hyper-angles associated with the directions of kx and ky. 
The three-body energy is then defined by 


The values Jijk depend on the volume in which the three- 
body systems are confined, and to obtain the a-*), coeffi- 


2 m 


( 27 ) 
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The arguments of the two hyperspherical harmonics, 3^^, 
are related to the spatial and momentum coordinates, 
respectively. The collection of angular quantum num¬ 
bers are denoted /C, where we only need to specify the 
hyper-momentum quantum number, K. The kinetic en¬ 
ergy eigenvalue corresponding to the free wave-function 
is K{K -I- 4) and at asymptotically large values of p 
related to the eigenvalue of the adiabatic potentials as 
K = i/{p —>■ oo) — 2. The asymptotic value of v for the 
lowest adiabatic potential is 2, corresponding to JtT = 0. 

The radial dependence is given in terms of the spher¬ 
ical Bessel function, Jk+ 2 , of order K + 2. The hyper¬ 
harmonics are normalized to unity 

J dny^,y,c = Sk'k, (28) 

where the delta-function express that all quantum num¬ 
bers pairwise must be equal. The large-distance asymp¬ 
totics, p — >• oo, is then obtained from the Bessel function, 
that is 


where we do not have to distinguish between choice of 
Jacobi coordinates, since the mass factor is symmetric 
under exchange of i, j and k. To take into account the 
degeneracy of the initial states for a given k we need to 
average over which we do by integrating over the 
angles and dividing by / dfln = iJs = tt^. For a given k 
(or energy) , and under the assumption that we are at 
large values of p, we then get 


J = 


/ my — 


■mj +mk^3/2dn^ 5 

—^-) -T^P dnpAjp. 34 

mimjmk 1^5 


Exploiting the orthonormality of the hyperspherical har¬ 
monics, the missing probability per unit time is then 


^pm, + mj+mk~^3/2 1 hn f 2TrY 

■' = “ ( > OT m (t) )i75- 

(35) 


In order to obtain the recombination coefficients we must 
multiply by V'^. This gives 


JK+ 2 iKp) 


p —^00 


2 ^IKp—iipK _j_ ^ — iKp+itpK 
TTKp 2 


(29) 


where ^pk = (^ + 2)f — f. The two terms in Ea. d^ 
correspond to in- and out-going hyper-sperical waves, re¬ 
spectively. We now introduce a short-range optical po¬ 
tential and a unit amplitude on the in-coming wave. The 
absolute value of the out-going amplitude is then asymp¬ 
totically allowed to differ from unity. From Eas. (E51) and 
(|2^ we then obtain the asymptotic from the K = Q wave 
function, that is 


1 (2^)^ 
V {kpY 


yo{n^)yo{np) 


2 

TTKp 2 


(30) 

where 0 < l^nl < 1, since the optical potential acts as 
a sink of probability. The hyper-radial current is defined 

by 




+mk^3/2^^2 (l-|‘S'iip) 


rriimjmk 




(36) 

which corresponds to the general formula for N-body loss 
given in [s^. 


C. Decay-rate method 

In this section we present an alternate way to obtain 
the recombination rate J for the three-body system. This 
is done by deriving the decay rate of bound states in 
the optical potential. First we define the bound states 
in a large box with hyper-radius extending from zero to 
Pmax- The boundary condition is that the wave function 
is zero at the edge of the box, that is f{pmax) = 0. The 
eigenvalues for the optical potential are complex numbers 


jp — 





(31) 


E — Eq — 



(37) 


Using Eas. (l30)) and (1^ the missing current, Ajp = 
jpiSii = 1) — jp{Sii), is calculated for a given value of 
S'!!. This amounts to 

(32) 

To get the probability loss per unit time, J, we need to 
integrate the missing current over the hyper-surface. The 
surface element must correspond to the physical volume 
element (see [H). The volume element is thus given by 
(see Eqs. dll) and ©): 


/ ^ N 3/2 

(fr^d^ry = (-] d^xd^y = 

\ PiPjk J 


2 ,,mi + mn + mk 

mi - - - 

mimjmk 


f'%^dpdnp, 


(33) 


The imaginary component of the energy describes the de¬ 
cay rate of the probability as seen from the time evolution 
of the wave-function, defined by 

\f{p^t)? = I/(pU = 0)|^exp(-rt/?i). 

The decay rates described by F = J7i are three-body 
energies determined by a box boundary condition. The 
overall energy dependence is then a strong decrease, in¬ 
versely proportional to the hyper-radial three-body vol¬ 
ume, towards zero as function of box radius pmax- In or¬ 
der to obtain the recombination coefficient f we then 
need to know V. This volume, V, is defined by equating 
two ways of calculating the density of three-body states 
in the hyper-spherical box extending to Pmax- The first is 
the formal expression of integration over given intervals 
of coordinates and conjugate momenta, where Px = ?ikx 
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and Py = Kky. The second is direct numerical calcula¬ 
tion of the same quantity from the solutions to the hyper- 
radial equation. The resulting equation is then 


J (fxd^yd^p^(fpyS{E - {pi +pl)/{2m)) = 

dn/dE = {2^%)^ dv / dE (38) 

where = f d^xd^y is the volume for the Jacobi co¬ 
ordinates, fl^K^dn/dE is the volume element per unit 
energy in momentum space, and du/dE is the density 
of states for a given energy, E, defined by Ea. (l27l) . The 
factor 2TTfi = h is Plancks constant, which is the volume 
occupied by each quantum state. 

Using Eg. (1551) we express V' in terms of the physical 
volume, y, that is 


y 2 = 


d'^r^d’^Ty = 


^,2^s^ m,+mj+mk y/2 _ 
^ miirijmk ^ 


Using Eas. (l551) . (1571) and (1551) we then get 

y 2 _ ^3 f "mi + mj + nik ^3/2 1 dv 2 /?i^(27r)^ 

'' rmmjmk E^ dE U 5 \ 2m 

(40) 

where the value of ^ can be found from solving the 
hyper-radial equation numerically. 

We emphasize that the lowest hyper-radial equation 
only accounts for both total and partial-wave angular 
momentum zero states, while the employed phase-space 
identity includes all states, independent of angular mo¬ 
mentum. The derived relations therefore strongly assume 
excitation energies sufficiently small to exclude contribu¬ 
tions from all solutions build on the repulsive higher-lying 
adiabatic potentials. 

The squared volume scales as oc 
meaningful recombination coefficients, independent of 
box size, are therefore only achieved when T is propor¬ 
tional to 2 ,. This is very demanding for numerical cal¬ 
culations, since convergence only is achieved when pmax 
is larger than the scattering lengths. 

The recombination coefficients, expressed in terms of 
r, are then found from 



tta 2 3.mi+mj +mk 

JV = m ( - - - 

mimjrrik 



1 dv ^ 
£"2 dE 


(41) 


Eg. dTTI) is equal to Eg. (1561) then. This leads to the im¬ 
mediate conclusion that for an energy E corresponding 
to an allowed eigenvalue the following relation between 
the parameters in the two methods should hold 


r = ^^(i-|^iip). (42) 


D. Recombination coefficients and finite 
temperature 


One immediate conclusion which is readily obtained 
from Eg. (1551) - Eg. (155|) is that 


,,( 2 ) 


_ 1 ( 1 ) 


^221 


*221 5 


(43) 


which is reassuring, as for example one particle of type 
2 disappears for each particle of type 1 in the 1-1-2 re¬ 
combination. For identical particles the mass factor in 
JU^ fEo. (156l) and Eo. (l4T1) l reduces to where rm is 
the physical mass. Numerically, we find that the tun¬ 
nelling probability 1 — 151112 (for an analytical estimate 
of this, see |^) and the decay rate T depends on 
the physical mass as mf, however, leading to the iden¬ 


tity a\^ii ~ ^® 222 - It i® approximate since the Efimov 
peaks can have different locations and shapes in the two 
systems, leading to non-systematic differences between 
r or equivalently 5ii in the two systems. This mass- 
dependence for identical bosons is in correspondence with 
earlier work [s^- 

In order to compare our calculations with experimental 
data we need to fold our calculated values of a[%{E) with 
a temperature distribution. The normalised Boltzmann 
distribution for 3 particles is given by [3l| 


= 2(W / E)dE , (44) 

where the E^ factor arises from the three-body phase- 
space. In order to get good results it is necessary to 
calculate in a range of energies around ki,T, where 
the integral receives contributions. 

The highest excitation energies allowed in our low- 
energy model are given by the energy difference between 
states build on the first and the neglected second and 
higher adiabatic potentials. This energy difference can 
be estimated by the difference in potential energies at 
distances where the states are located. Since the scatter¬ 
ing length is a measure of the sizes of all these Efimov- 
like states, we use the hyper-radius p ~ IO^oq to give a 
lower estimate of the maximum temperature allowed in 
realistic calculations. The energy difference between first 
and second adibatic potential is for interaction free states 
given by fi^K{K + A)/ {2mp^) where K = 2 The re¬ 
sult of these estimates are temperatures of the order of 
~ IpK, which therefore is an upper temperature limit 
for realistic calculations. 

With the temperature distribution implemented by 
Eg. (1551) the two methods, that is the S-matrix and the the 
decay rate, give essentially the same results for the same 
choice of pcut and Vimag, aside from numerical inaccura¬ 
cies. We shall use whichever method is the most conve¬ 
nient in the practical calculations. For technical reasons 
it is for example generally more convenient to implement 
the zero-energy limit by the decay rate method, while it 
is easier to implement the temperature distribution for 
the S-matrix method. 







IV. PARAMETER DEPENDENCE 


0.5 


The recombination coefficients for our two-component 
system depend on three parameters, that is one mass 
ratio and two scattering lengths. We use the notation 
from fig. [T] where the particles 1 and 3 are identical, 
and 2 is distinctly different. To illustrate the effects we 
investigate first the variation of the adiabatic potential 
which is the crucial ingredient in all the calculations. Be¬ 
fore we investigate the dependence of the recombination 
on physical parameters, we investigate the dependence 
on the optical model parameters. Then we investigate 
the dependence on scattering lengths for masses of sys¬ 
tems where experimental results are available. Finally 
we compare the different recombinations that can occur 
in a two-component gas. The factors corresponding to 
a non-BEC gas are used for all the recombination coeffi¬ 
cients. 

The numerical calculations are technically simple and 
employ only homemade standard programs. First the 
lowest adiabatic potential is found numerically from the 
complex angular eigenvalues, which in tur is found by 
solving Eq.®. The absorption probability is then calcu¬ 
lated for different energies from the probability reduction 
of a plane wave reflected by the optical model poten¬ 
tial. This involves solving Eq.(IS]), at a given energy, for 
the modified potential Eg. (11^1) . This is done by numeri¬ 
cal integration with initial conditions in pcut provided by 
Eg. lfTSl) . The result are then fitted to Ea. (l30l) in order to 
extract i'll. The complex eigenvalue of the optical po¬ 
tential is calculated by the shooting method, by requiring 
that the wave-function is zero at Pmax- Tfie numerical 
integration is implemented as in the above i-matrix cal¬ 
culations. 


A. Adiabatic potentials 

The masses only enter as ratios of masses through the 
(f>ij functions in Eq.®, that is for our case this leaves 
only one parameter, R = ^. We show computed mass 
dependence of adiabatic potentials in fig.[3]as functions of 
hyper-radius measured in units of the scattering length, 
ai 2 of the distinguishable particles. The two figures show 
results for a vanishing and a relatively large scattering 
length between the identical particles an. We see that 
all potentials asymptotically approach , as in the case 
of identical particles. 

For both values of an we find that increasing R leads 
to lower barrier height, and a shift to smaller hyper-radii 
of barrier position and p-value, where the potential is 
zero. The dependence is strongest for two heavy and one 
light particle, that is for R > 1. Continuing to i? < 1 
would give very little variation in the potentials compared 
to the R = 1 curve. Increasing |an| from zero reduces 
the peak heights and move the potentials towards larger 
hyper-radii. 



FIG. 3: The lowest adiabatic potentials as function of p/ai 2 , 
calculated from Eqs. m and m- The lower figure is for 
an = 0 and the upper figure is for on = 012 . 



FIG. 4: The lowest adiabatic potentials as function of 
p/oi 2 . The red graph corresponds to Ea. (ll2ll . The 
other graphs, in descending order, corresponds to an = 
0.08, 0.16, 0.32, 0.64,1.28ai2 where 012 is negative, calculated 
from Eq. Hill- 


Experimental data are available for a Cs-Li gas where 
the mass ratio is R = 22.28. The dependence on the 
scattering length, ai 2 , is shown in fig. |T] for different val¬ 
ues of an. The potential is rather small for this value 
of R, but still with a distinct barrier where the height 
decreases with increasing values of |aii|. 

To complement we show the dependence on an in fig. [5] 
for different values of 012 . The barrier decreases with 
increasing 012/011 from a rather large value for vanishing 
|ai 2 |, towards rather small values. For large values of 
|ai 2 | the barrier has become so small and moved so far 
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FIG. 5: The lowest adiabatic potentials as function of pjaw. 
These graphs, in descending order, corresponds to ai 2 = 
0.01, 0.02, 0.04, 0.08, 0.16, 0.32, 0.64, I .28011 where an is neg¬ 
ative, calculated from Eq. (El). 

to the right, that it is no longer visible on the scale chosen 
for the figure. 

These figures show that increasing |aii| and |ai 2 | has 
the same qualitative effect as increasing the scattering 
length in the case of identical bosons for the optical 
model, see fig. [2] It decreases the height of the barrier 
and moves the location of its maximum to the right. 

B. Optical model 

The real potential is completely determined for zero- 
range two-body interactions in terms of scattering 
lengths and masses. The box-like imaginary potential 
has hyper-radius and strength as the two phenomenolog¬ 
ical parameters. This radius can only assume discrete 
values determined by the requirement that the scatter¬ 
ing length is reproduced for the occurrence of a specific 
recombination peak. The corresponding strength corre- 
latedly varies shape and size of the chosen peak. On fig.jS] 
the results of varying Vimag, while pcut is a constant value 
of 1.4ao are shown. 



FIG. 6: The recombination coefficient for Cs-Gs-Cs as func¬ 
tion of the scattering length a for different choices of Vimag, 
at pcut 1.4ao 

We see that increasing the size of Vimag makes the 


peaks more pronounced while also making the absolute 
value for the rest of the graph somewhat smaller. De¬ 
creasing the size of Vimag has the opposite effect. This 
is the exact same qualitative effect which was seen when 
varying the absorption parameter n* describing deeply 
bound states in 

From one set of parameters another equally good 
choice is found by scaling the square well hyper-radius 
up or down by the Efimov factor and at the same time 
scaling the strength in opposite direction by the square of 
the same Efimov factor. Using these parameters repro¬ 
duces exactly the same shape of the Efimov resonance 
and the same absolute value of the recombination. This 
means that Vimag pint fo'' ^ given peak has a constant 
value where pcut can discretely vary by the Efimov fac¬ 
tor. The largest allowed hyper-radius must be at a hyper¬ 
radius where the real potential precisely has the inverse 
radial square behavior. The lowest value is only limited 
by the requirement that it must be finite, since otherwise 
the already removed zero-range divergency reappears. In 
practical calculations we generally pick values in the or¬ 
der of oo to be certain that we do not approach any of 
these limiting cases. 


C. Recombination 

For identical particles we know that the recombination 
coefficient is proportional to the fourth power of the scat¬ 
tering length. For a two-component system this relation 
must be extended to account for two different scattering 
lengths as well as for a mass-ratio dependence. We shall 
illustrate with two mass ratios R = 22.2 and R = 2.23, 
corresponding to experimentally realizable systems, that 
is ®Li-^^^Cs 0, im and ^®K-®’’Rb gases. We used these 
isotopes to have well-defined mass ratios in the calcula¬ 
tions. The characteristic small mass variations between 
isotopes would only marginally change the results, pro¬ 
vided the boson-fermion characters remain the same or 
become irrelevant as when only one identical atom par¬ 
ticipates in the process. Here we shall only investigate 
the pure mass dependence. 

The dominating recombination coefficient is related to 
= 20 ^ 22 j where label 2 corresponds to the light par¬ 
ticle. This means that we here consider the heavy-heavy- 
light three-body systems with mass ratios, R = 22.2 and 
R = 2.23. We calculate all the recombination coefficients 
in the limit of zero three-body energy. 

The periodic structure of enhanced recombination oc¬ 
curs each time the scattering length is multiplied by the 
Efimov scaling factor, s = exp(7r/|i/(/5 = 00 )|), found for 
infinite scattering lengths of the contributing systems. 
This scaling parameter depends first of all on the mass 
ratio. The Efimov effect requires that |ai 2 | = 00 , while 
|aii| can assume any finite or infinite values. The results 
for the two limiting cases, |aii| = 0 and 00 , are given in 
tableUas a function of mass ratio, R. For large R the two 
cases are almost identical, but there is a big difference for 
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small mass-ratios and this trend continues for R < 1. If 
all three scattering lengths are infinitely large the scaling 
for i? = 1 is s = 22.7 and for small R approaches a con¬ 
stant of s = 15.7. For a more detailed discussion of the 
two different cases see [2^. 

TABLE I: Efimov scaling factor s for different values of R 


-R = ^ 

m 2 

1 

2 

5 

10 

15 

20 

ai 2 = dll = 00 

22.7 

20.8 

13.7 

8.50 

6.30 

5.14 

ai 2 = 00 , ail =0 

1986 

153.8 

23.3 

9.76 

6.64 

5.25 


The influence of an can conveniently be studied 
through the recombination of the two chosen systems. 
For R = 22.2 the Efimov scaling only varies between 
4.8766 for an = 0 and 4.7989 for |an| = oo and the peak 
positions should therefore remain. In contrast the scaling 
between peaks for R = 2.23 should move between the ex¬ 
treme limits of 121.1 for an = 0 and 20.28 for |an| = oo 
as a function of an. 



FIG. 7: The recombination coefficient, an 2 > function of 
ai 2 , for different values of an, where R = 22.2. The param- 
eters of the optical potential are Vimag = with pent 

in the interval, 0.24-0.32ao as an varies from 0 to -20000ao, 
adjusted to maintain the position of the first Efimov peak at 
ai 2 ss —320ao. 

Let us first focus on the recombination coefficient, ctnii 
for a large mass ratio of 22.28. We calculate 0^2 

a function of ai 2 , for different constant values of an. 

*2 

The strength of the optical potential is Vimag = 68^^, 
whereas peut is adjusted slightly to reproduce the peak 
position |aj Ri — 320ao, corresponding to the experi¬ 
mental peak position of Cs-Cs-Li [20,[^, for the different 
values of an. 

The main results of these calculations are summed up 
in fig. [71 where we show the coefficient as a function of 
ai 2 for different values of an. The most striking fea¬ 
ture of fig. [7] is that the different values of an lead to 
a different overall dependence of a[^i 2 as function of ai 2 . 
For small values of |an| there is a a^ dependence cor¬ 
responding to the scaling for identical bosons, but for 
bigger values of |an| this relation is no longer valid. As 
|an| grows bigger, the recombination is enhanced, which 


is most visible at smaller values of |ai 2 |. This is consis¬ 
tent with larger values of each of the scattering lengths, 
|ai 2 | and |ai 2 |, lowering the potential barrier, making it 
easier to reach small distances and thereby enhancing the 
recombination. 

We also know that ai 2 has a stronger influence on the 
potential than an, which corresponds well with ai 2 be¬ 
ing the most important parameter for the recombination 
shown in fig. [71 Specifically, we see that 0^2 changes 
more as ai 2 runs from 0 to -20000 ag, than it does when 
an varies from 0 to -20000 uq. 

All in all, for realistic values of the scattering lengths, 
we can view an as moderating the ai 2 dependence and 
simultaneously enhancing the recombination. More vio¬ 
lent changes of an cannot be excluded when the interac¬ 
tions are controlled by the Feshbach resonance technique. 
However, this is not the case in any of the physical sys¬ 
tems discussed here. 



FIG. 8: The recombination coefficient, Ojjj s-s a function of 
ai 2 calculated for different values of an, where R — 2.23. The 
strength of the potential is as in fig. [7] while pc.ut is adjusted 
to give one peak at about ai 2 « —300ao. 

We now move to the mass ratio R = 2.23 where the 
Efimov scaling has a strong dependence on an as seen 
from the extreme limits given in table [H The relative 
positions of the recombination peaks must then vary with 
the finite value of an. Again we calculate an 2 ^ 
function of ai 2 for different values of an. The results of 
such a series of calculations are shown in fig. [51 where 
Pent is adjusted to produce the same somewhat arbitrary 
peak position at about |a^ r; —300ao for all an. The 
next Efimov peak then has to move in the interval of ai 2 
from 20.28-a^^i R —6000ao and 121.1-ai^i R —36000ao. 

TABLE If: The ratio, \ between the ai 2 scattering 

lengths for the first two peaks as function of an when R — 
2.23. The last column give the related peut values. 


an [ao] 

0 

-1000 

-5000 

-10000 

-20000 

-50000 

TFT 

119.4 

96.9 

41.9 

32.9 

27.7 

23.9 

Peut [o-o] 

1.21 

0.97 

1.14 

1.18 

1.2 

1.22 


The most noticeable thing in fig. |5| is that the loca- 








































11 


tion of the second Efimov resonance changes with the 
finite value of an. This variation of the ratio between 
ai 2 values of second and first peak is shown in table HIl 
for different an values. The rather modest variation of 
Pent is also shown in this table. We notice the correct 
continuous variation between the two extreme limits of 
the Efimov scaling with an, although the value of peut for 
an = 0 reflects the singularity moving between Eqs. (nni) 
and IT^ . 

We emphasize that small changes of Peut is fine-tuning 
to maintain the first peak at the same position. The 
much larger shift of the second peak is entirely due to 
the variation of an, and essentially completely indepen¬ 
dent of the parameters of the imaginary potential. A 
substantial change in the second peak position requires 
a value of an ~ —5000ao. However, the change is fairly 
gradual and there is no ’’magic value” where an sud¬ 
denly begins to contribute. The ideal Efimov scaling for 
three resonant interactions of 20.28 is essentially reached 
numerically when an = —50000ao. 

The second prominent feature in hg. [8] is that a finite 
value of an modifies the overall dependence of ^ 

function of ai 2 . This is the same qualitative feature as 
seen in fig.[7]for the large mass ratio R = 22.2. Einite val¬ 
ues of Oil leads to a higher absolute value of cx^n 2 within a 
huge ai 2 interval, ai 2 G [0, —20000ao]. This means that 
increasing the size of |aii| enhances the recombination 
coefficient, like we found for the R = 22.2 system. 

D. Different recombination processes 

We have so far only considered recombination from 
three-body systems with two heavy and one light parti¬ 
cle such as Rb-Rb-K and Cs-Cs-Li. However, the same 
two-component gas can also decay by the other 3 combi¬ 
nations, two light and one heavy particle (Rb-K-K, Cs- 
Li-Li), and 3 identical heavy (Rb-Rb-Rb, Cs-Cs-Cs) or 
light particles (K-K-K, Li-Li-Li). In order to estimate the 
importance of the terms in Eas. (fl6l) and ini) , it is nec¬ 
essary to calculate all these recombination coefficients. 

To make a realistic comparison we first consider the Cs- 
Li gas, where the Efimov resonances are fixed by experi¬ 
mental data for Cs-Cs-Li [2^, Cs-Cs-Cs [l^ and Li-Li-Li 
[l^ . and the parameters of Li-Li-Cs are assumed to be 
the same as for Cs-Cs-Li. We assume that an = 022 = 0. 
This would underestimate the recombination coefficients 
as shown in the previous subsection, but it allows on the 
other hand a clean comparison where the fourth power 
dependence applies, oc a^- Actual correlated finite 
values of these scattering lengths obtained through the 
Feshbach technique could then be important and quan¬ 
titatively alter the comparison. 

On hgure [5] all the recombination coefficients for the 
Cs-Li gas have been plotted. The optical model param¬ 
eters are chosen to reproduce the position of the lowest 
measured recombination peaks [171 Il8l . l20ll . The exper¬ 
imental data for Li is actually for ^Li |l7j |. but our aim 


is here only to test the mass dependency. In the ^^^Cs- 
®Li experiment, Li-Li-Li and Li-Li-Cs recombinations are 
expected to be suppressed due to Fermi statistics, which 
is not taken into account in our model. So the results 
are not directly comparable to the experiment. In ad¬ 
dition this is not a direct comparison between the mass- 
balanced and mass-imbalanced case, since the recombina¬ 
tion coefficient is a function of different scattering lengths 
in the different cases. 

Furthermore, for 022 = 0 it was not possible to lo¬ 
cate an Efimov resonance within the range of ai 2 € 
[0, — lOOOOoo] for neither the Li-Li-Cs nor the K-K-Rb 
mass ratio, even though a wide range of peut were tested. 
This is because the Efimov scaling factor now is so big 
that it is hard to locate an interval with even one Efi¬ 
mov resonance. In these comparisons we shall focus on 
recombination coefficients for zero energy. 



FIG. 9: The top shows the recombination coefficients result¬ 
ing in Cs loss, and the bottom shows the recombination co¬ 
efficients resulting in Li loss, both as a function of ai 2 . The 
optical parameters in the calculations are peut = 0.20ao and 
Vimaq ~ 240-L-- is used for Cs-Cs-Li and Li-Li-Cs. For the 

^ mag 

t 2 

Cs-Cs-Cs recombination pent = 1.4ao and Vimag = 16;^^ are 

used and for Li-Li-Li pent = 0.41oo and Vimag = 68^^ are 
used. 

The absolute sizes on fig. [9] show that the Cs-Cs-Li 
recombination is much more likely than the Li-Li-Cs re¬ 
combination for the same value of ai 2 . The Cs-Cs-Cs 
recombination only depends on an, which is used as the 
x-coordinate for this process in fig. HI This scattering 
length is expected to be of less importance compared to 
ai 2 in the mixed recombination coefficients, and therefore 























12 


assumed to be zero in those estimates. 

The comparison is then not straightforward but still 
useful, since a finite value of an would increase aJn b®" 
yond the curve in fig. HI Thus we deduce that aj^ 2 (®ii = 
0 , 012 ) ^ = “ 12 )) and we therefore believe that 

the Cs-Cs-Li recombination is much more likely to oc¬ 
cur than the Cs-Cs-Cs recombination in realistic systems. 
We also see in fig. H] that q.\ 22 { 0'22 = 012 ) ^ = 

0, 012 ). This does not allow any conjecture about relative 
sizes in a realistic system because the intra-species scat¬ 
tering lengths usually are much smaller than the neces¬ 
sary large (for the Efimov effect) inter-species scattering 
length. The recombination coefficients between identical 
particles are then expected to be relatively small. 



FIG. 10: Recombination coefficients for R = 2.23 as a func¬ 
tion of a 12 . The optical parameters are pent = 2.79ao and 
12 

Vimag = 20^^ for both systems. 


cated relations between the different recombination coef¬ 
ficients which only can be determined by taking the Fes- 
hbach resonances of a specific system into account. So in 
theory, it is possible to have a specific system where the 
light-light scattering length |a 22 | is much bigger than the 
heavy-heavy scattering length |aii|, which in turn leads 
to dominance of the light-light-heavy over the heavy- 
heavy-light process. 


V. COMPARISON WITH EXPERIMENT 

In this section we will confront our theoretical results 
with recent experiments that have been carried out for a 
^^^Cs-®Li gas [20,[2l|, as well as with a recent experiment 
for a ^^^Cs gas, in which a second Efimov resonance has 
been observed The variable parameter in experi¬ 

ments is the magnetic field B, which via the mechanism 
of Feshbach resonances can be used to change the scatter¬ 
ing lengths, a. The phenomenological relation between a 
and B is 


a{B) = Qbg (^1 -b 

where A , abg and Bq are determined experimentally for 
each individual system. Eg. (1451) applies for most systems 
and will be used in the cases investigated in this section. 
We first focus on the equal mass process after which we 
continue with the dominating heavy-heavy-light recom¬ 
bination process. 


To complement we now investigate the much smaller 
mass ratio R = 2.23 (corresponding to a Rb-K gas). The 
results are shown on fig. (TU] for different recombination 
processes. The optical model parameters are chosen to 
give a peak for ai 2 ~ — 700ao in the coefficient. 
As for the larger mass ratio we again conclude that the 
heavy-heavy-light recombination process is more likely 
than the ligt-light-heavy recombination process. 

Figs. H] and [TO] also lead to another conclusion. A 
bigger mass-ratio seems to give a bigger recombination 
coefficient for two heavy particles and one light. This 
is in accordance with a corresponding decrease of size of 
the adiabatic potentials, which intuitively suggests higher 
probability for recombination at small hyper-radii. 

A bigger mass-ratio also seems to give a somewhat 
smaller recombination coefficient for two light particles 
and one heavy particle, although not as pronounced. 
This means that the bigger the mass-ratio, the bigger the 
difference between the heavy-heavy-light recombination 
and the light-light-heavy recombination. For both small 
and big mass-ratios, the 0:221 terms can be neglected in 
Eqs. (O and (II 3 based solely on these mass-related ar¬ 
guments. 

We emphasize these conclusions are only strictly valid 
for small jan | and |a 22 |. Finite values imply more compli- 


A. The recombination 

Recently a second Efimov resonance was observed in a 
133 Cg gas and we can conveniently test the model 

against this experimental confirmation of the original Efi¬ 
mov scenario. The first Efimov resonance for a ^^^Cs 
gas was observed earlier [l^ to have a peak for Ri 
—960ao. The experimental data from these measure¬ 
ments also allow us to adjust both peut to give the peak 
position and subsequently tune the strength, Vimag, to 
the shape of this first peak. With the scaling mass, m, of 
^^^Cs, peut = l-4ao and Vimag = 16;^, the first Efimov 
resonance is then rather well reproduced. 

We show the calculated results on fig. [11] where the 
second Efimov peak is obtained without any adjustments 
beyond the first peak. The experimental results are not 
completely commensurable, since the first experiment 
was done at a temperature of about 15 nK, and the 
second experiment at 9.6 nK. We can circumvent this 
by using the fact that temperatures of this value are in¬ 
significant at small values of the scattering length. The 
data for the first peak from the first experiment [l^ can 
therefore be assumed to arise for the same temperature 
of 9.6 nK as the second peak in the second experiment 
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FIG. 11: The recombination coefficient for Cs-Cs-Cs, ajVn 
different temperatures plotted with experimental data from 
Thfi optical model parameters are pcut = 1.4ao and 
V- — 16 


0- Both temperatures are much smaller than the up¬ 
per temperature limit estimated to be IfiK for a realistic 
calculation. 

The absolute value around the first peak and for 
big a is generally remarkably close to the experimental 
value. However, the calculated recombination coefficient 
is much too small at small a, where the temperature has 
no influence. The theoretical results follow the overall 
rule of recombination dependence and these small a 
deviations must therefore arise from other processes con¬ 
tributing to the experimental values. 

The calculated temperature dependence for larger 
scattering lengths is in very good agreement with the 
measurements. This is remarkable, since a reduction in 
temperature to 1 nK results in a fairly dramatic change 
of the recombination curve at large values of a. The con¬ 
sequence is that a moderate temperature of a few nK al¬ 
ready smears out the second Efimov peak, and prohibits 
observation. We then conclude that using the correct 
temperature gives a shape that is in pretty good agree¬ 
ment with the experimental results. However, the second 
peak is dislocated compared to the ideal Efimov scaling 
factor of 22.7, which predicts this peak to be at around 
ai 2 = —21800 oq. In the experiment, it is found at ap¬ 
proximately — 17000 oo, while the calculation gives the 
peak at —15800 oq. 

Overall, the results of this comparison with experi¬ 
ments for identical bosons are encouraging. The tem¬ 
perature effects are well accounted for and the shape and 
location of the second Efimov peak is also in broad agree¬ 
ment with data, for phenomenological parameters fitted 
to the first peak. Our model seems to work for the well- 
known mass-balanced case, and comparison to data for 
mass-imbalanced systems should then be considered. 


B. The ^^^Cs-^^^Cs-^Li recombination 

The crucial parameter is the scattering length, ai 2 = 
flLzCsj which has to be very large to provide the Efimov 
effect. However, in addition also an = acsCs is impor¬ 
tant for quantitative predictions. The overall scaling 
is modified for a finite value of an, which is determined 
by the magnetic field through the Feshbach resonance of 
the system as described in Ea. (H5|) . It is then interest¬ 
ing to look for effects in the recently obtained two sets 
of experimental data HH . They are in broad agree¬ 
ment, although the details are a little different. In [2l| 
3 Efimov resonances are reported, where the third one 
is very hard to distinguish from the background due to 
finite temperature effects. In the recombination coef¬ 
ficient is given as a function of ai 2 , which allows an easy 
comparison with our calculations. 

The experimental conditions in provide the pa¬ 
rameters in Ea. (H51) for the interspecies (Cs-Li) Feshbach- 
resonance of the prepared spin-states, that is 

Qbg = -28.5ao , A = 61.4 G , Hq = 842.9 G , (46) 

which by insertion in Ea. (l45|) gives the variation of ai 2 . 
The Oil scattering len gth is estimated to run between 
— 1200ao and — 1500ao so we have chosen a constant 
intermediate value of an = acsCs = —1350ao. 



FIG. 12: The recombination coefficient for Cs-Gs-Li,Q ()2 as a 
function of the interspecies scattering length, for an = 0 and 
an = — 1350ao, plotted with the experimental data from [^ . 
The optical model parameters are pcut = 0.258ao , Vimag = 

fc2 

100^^ for an = -1350ao and pcut = 0.20ao , Vimag = 
240-i^7 for an = 0. One small temperature of 1 nK and the 

ma-Q ^ 

experimental temperature of 450 nK are shown. 

We compare theoretical and experimental results in 
fig. [T2] where both temperature and an dependence are 
shown. We hrst note that there is a pretty good corre¬ 
spondence between the experimental data at 450 nK and 
the theoretical prediction at 450 nK with the suggested 
finite value of an- The shape is nicely reproduced and 
the absolute value is within the systematic uncertainty. 
We emphasize that the calculations for vanishing an = 0 
does not describe the data quantitatively at 450 nK. The 
absolute value and the shape is only correct when the 
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enhancement of a^ii 2 at small values of ai 2 , due to finite 
values of an, is taken into account. We also draw atten¬ 
tion to the fact that for a substantial finite value of |aii|, 
Q;n 2 longer independent of energy at small values 

of |ai 2 |. The temperature of 450 nK is not far from l^K, 
that is upper temperature limit in a realistic calculation. 
But the data is described very well. 


It seems like the calculated values of describe the 
relative behaviour of different recombinations well. The 
recombination coefficient for Cs-Cs-Li is larger than for 
Cs-Cs-Cs for both the experimental and the calculated 
values. The calculated values of both are within experi¬ 
mental uncertainty. The relative location of the Efimov 
resonances depends on the finite value of an. For the 
Cs-Cs-Li system the limiting values are 4.8766 or 4.7989, 
however, which means that we do not expect any note- 
worthy difference. The 3 Efimov resonances in 21| are 
well suited for testing the theoretical prediction for the 
Efimov resonances. 

We use the value, Bq = 942.75, of the resonant mag¬ 
netic field. The two reported different values of Bq are 
within the experimental uncertainties of both experi¬ 
ments. In table uni we compare predicted and exper¬ 
imental values for the Efimov peaks when the optical 
model parameters are adjusted to reproduce either first 
(pcut = 0.258ao) or second peak {pcut = 0.28ao). 


TABLE III: The values of AB — B — Bq determined to give 
the first or second Efimov peak for Cs-Cs-Li, while the other 
two peaks ( second or first as well as third) are predicted with 
that set of corresponding optical model parameters. 



1st peak 

2nd peak 

3rd peak 

Experimental value of AB [G] 

5.61(16) 

1.07(2) 

0.22(4) 

AB [G] fitted to 2nd peak 

5.8 

1.07 

0.21 

AB [G] fitted to 1st peak 

5.65 

1.03 

0.2 


Then in both cases the third peak is in agreement with 
calculations within the experimental uncertainty. How¬ 
ever, this does not say much, since a small change in 
magnetic field results in a huge change in the scatter¬ 
ing length, see Ea. (H51) and Ea. (H51) . Due to this, it is 
also hard to determine an accurate pcut value which re¬ 
produces the very weak third peak. The two predicted 
peaks in the two fits are at the edge allowed by the ex¬ 
perimental uncertainty, and as such they are still very 
close to the observed values. The finite temperature ef¬ 
fects are also likely to move the locations slightly, which 
might explain smaller deviation. 

We can now calculate the inter-species scattering 
lengths corresponding to the experimental peak values 
by using Ea. (l46)) with Bq = 942.75, which yields |a^~^| = 

331. 5ao, |a 2 = I62Iao and jag = 7777. 5ao. This 

gives a Efimov scaling of 4.8911 between the first two 
peaks and 4.7980 between the last two peaks. These val¬ 
ues are in agreement with the Efimov scaling found by 
going to the universal limits of infinite scattering lengths. 


4.8766 or 4.7989, depending on whether there is two or 
three contributing resonant interactions. 


C. Universal properties of optical parameters 


For identical bosons ref.[3l| proposes a possible rela¬ 
tion between the optical model parameters and the van 
der Waals length of the system. This is motivated by the 
recent findings in equal mass systems [l^ of a univer¬ 
sal relation between the three-body parameter and the 
two-body van der Waals length 4TI-I4^. We can expand 
this idea to the mass-imbalanced system. The three van 
der Waals lengths which we will compare with is for Cs- 
Cs. Li-Li and Cs-Li respectively 202ao, 65ao and 44.8ao 
[ 43 , The crucial parameter is pcut which by definition 
is a hyper-radius. From the definition in Eq.(IT]) we can 
analogously define a hyper-radial van der Waals length 
as 


2R 


Mil 


1 + R 


PvdW i^2R m "’"‘^^41 + 1/2 + i? m 


M12 2 

^vdW,l2 


(47) 


where the distance between particles is replaced by the 
corresponding two-body van der Waals lengths, and the 
related two-body reduced masses, pn and pi 2 , are intro¬ 
duced and R= 

For identical particles with m equal to the mass of the 
particles, this hyper-radius reduces to the van der Waals 
length, which in ref.[3l| was compared 

to Pcut by simple division. For three identical particles we 
find respectively Pcut/Pvdw = 0.0069 and Pcut/Pvdw = 
0.0063 with the parameters for Cs-Cs-Cs and Li-Li-Li. 

For the mass-imbalanced system the (identical) inter¬ 
actions between the heavy and the two light particles is 
as dicussed expected to be the dominating contribution. 
This means that the last term in Ea. dTTl) should be the 
largest while the first term should vanish when the cor¬ 
responding scattering length, an, between the identical 
heavy particles is relatively small. Furthermore, when 
an is finite and begins to contribute the recombination 
peak moves. To keep the peak location independent of 
an we therefore varied pcut slightly with an. This ad¬ 
dition to Pcut has to come from the first term of Eq. (iTTl) 
which should vary from zero to a given hnite contribution 
when I an I = 00 . Thus we can parametrize by 


PvdW 


2R pii |an| 2 

I + 2i?^ (/nao + |an|i?2/(^))'’’"'^’“ 


I -|- i? pi2 2 


(48) 


where we suggest to use /n = 10^. For equal masses and 
an —>■ 00 this reduces to Eq|47]for equal masses. For un¬ 
equal masses the R^ factor ensures that for an — 00 , the 
term still doesn’t become dominant, in correspondence 
with the moderate change in pc«t, even as big values of 
|an| were used numerically. We use f{R) = I, but per¬ 
haps a more complicated function is required. For the 
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experimental value, an = — 1350ao {pcut = 0.258ao), we 
then get Pcut/Pvdw = 0.0057 close to the experimental 
values for identical particles. 

The strength, Vimag, can also be compared to the 
natural unit for a van der Waals interaction, V^dw = 
which amounts to VimagIVvdw = 6.53 • 
10^,2.87 • 10^,2.03 • 10^ for Cs-Cs-Cs, Li-Li-Li, and Cs- 
Cs-Li, respectively. The value of ^Vimagplut is constant 
for a given location and shape in a specific system. The 
numerical values from our calculation are for this quan¬ 
tity 31.3600, 11.4308, 6.6564 for the systems Cs-Cs-Cs, 
Li-Li-Li and Cs-Cs-Li. 

The upshot of these comparisons is that a possible ap¬ 
proximate value of Pcut for a system can be obtained 
by Pcut ~ Q-QQQpvdw- In addition the value of Vimag 
to within a factor of two seems obtainable by Vimag ~ 
3 • Vf'Vudw- Finally the value ^VimagPcut in different 
systems are within the same order of magnitude. 

VI. DISCUSSION AND OUTLOOK 

We have developed a method for calculating the three- 
body recombination coefficient of mass-imbalanced three- 
body systems at negative values of the scattering lengths 
and for small energies. In order to do this we employed 
the hyper-spherical method, zero-range potentials and 
the Faddev decomposition to calculate the lowest-order 
adiabatic potential. 

In addition we formulated and explored the optical 
model to calculate recombination processes into deep 
dimers. This method introduces two phenomenological 
optical parameters, strength and range, which respec¬ 
tively determine location and shape of the Efimov res¬ 
onances. We then developed two methods for finding 
the recombination coefficients of both mass-balanced and 
mass-imbalanced systems. By means of the traditional 
S-matrix method and by a new method, where the de¬ 
cay rate of bound states in a box due to the presence of 
the optical potential is calculated. In general the results 
obtained with the two methods are essentially indistin¬ 
guishable, although differing in numerical inaccuracies. 
As they tend to complement each other, we choose the 
most convenient method in actual calculations. 

The model was tested against the experimental data 
on recombination coefficients in Cs-Cs-Cs and Cs-Cs-Li 
systems as functions of the Feshbach tuned scattering 
lengths. In both cases, after fitting the range and the 


strength of the imaginary potential to the first peak, the 
model was able to describe quantitatively the whole curve 
including the temperature effects. 

The two-parameter fits of strength and hyper-radius 
of the imaginary part of the optical potential are very 
efficient for both position and shape of the recombina¬ 
tion peaks. It is also remarkable that the absolute values 
of the experimental recombination coefficients are repro¬ 
duced within the experimental uncertainty. 

Using the developed methods we have reached a num¬ 
ber of conclusions. The main conclusions are that re¬ 
combination is dominated by the heavy-heavy-light pro¬ 
cess and mainly determined by the heavy-light scattering 
length. But the heavy-heavy scattering length is impor¬ 
tant for obtaining the correct value of the recombination 
coefficients, as it enhances the recombination probability 
and modifies the behaviour as a function of the heavy- 
light scattering length. For a large mass ratio the Efimov 
scaling is determined entirely by the heavy-light scatter¬ 
ing length, but for a smaller mass ratio finite values of the 
heavy-heavy scattering length becomes important. The 
Efimov scaling then moves continuously between values 
from two and three resonating subsystems. 

The focus of this paper has been the case where all 
scattering lengths are negative, but it could be just as rel¬ 
evant to investigate two negative, one positive or two pos¬ 
itive, one negative scattering length cases. Since a finite 
value of the heavy-heavy scattering length can substan¬ 
tially alter the absolute value of a given two-component 
recombination, it may be interesting to test the three- 
body recombination for the same system at different Fes¬ 
hbach resonances experimentally, in order to further con¬ 
firm this prediction. 

For small mass-ratios we find a big difference in the 
Efimov scaling between peaks as the heavy-heavy scat¬ 
tering length is increased. It may be fruitful to investi¬ 
gate the intermediate area between two and three reso¬ 
nant subsystems in more theoretical detail. In addition 
a Feshbach resonance in a two-component gas with small 
mass-ratio, in which the heavy-heavy scattering length is 
big in the same area as the inter-species scattering length 
would allow investigation of this effect. 
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